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ABSTRACT 

If the gas in filaments and halos shares the same velocity field than the lu- 
minous matter, it will generate measurable temperature anisotropies due to the 
Kinematic Sunyaev-Zeldovich effect. We compute the distribution function of the 
KSZ signal produced by a typical filament and show it is highly non-gaussian. 
The combined contribution of the Thermal and Kinematic SZ effects of a filament 
of size L ~ 5Mpc and electron density Uf. ~ 10^m~^ could explain the cold spots 
of 5 ~ — 200/iK on scales of 30' found in the Corona Borealis Supercluster by the 
VSA experiment. PLANCK, with its large resolution and frequency coverage, 
could provide the first evidence of the existence of filaments in this region. The 
KSZ contribution of the network of filaments and halo structures to the radiation 
power spectrum peaks around / ~ 400, a scale very different from that of clusters 
of galaxies, with a maximum amplitude l{l + l)Ci/27r ~ 10 — 25{fiKy, depend- 
ing on model parameters, i.e., as and the Jeans length. About 80% of the signal 
comes from filaments with redshift z < 0.1. Adding this component to the intrin- 
sic Cosmic Microwave Background temperature anisotropies of the concordance 
model improves the fit to WMAP Syr data by Ax^ — 1. The improvement is not 
statistically significant but a more systematic study could demonstrate that gas 
could significantly contribute to the anisotropies measured by WMAP. 

Subject headings: Cosmic Microwave Background. Cosmology: theory. Cosmol- 
ogy: observations. 
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1. Introduction 

At low redshift, about 50 % of all baryons have not been identified (Fukugita & Pee- 
bles, 2004). The highly ionized intergalactic gas, that has evolved from the initial density 
perturbations into a complex network of mildly non-linear structures in the redshift interval 
< 2; < 6 (Rauch 1998; Stocke, ShuU & Penton 2004), could contain most of the baryons 
in the Universe (Rauch et al. 1997, Schaye 2001, Richter et al 2005). With cosmic evolu- 
tion, the fraction of baryons in these structures decreases as more matter is concentrated 
within compact virialized objects. The Lya forest absorbers at low redshifts are filaments 
(with low HI column densities) containing about 30% of all baryons (Stocke et al. 2004). 
Hydrodynamical simulations predict that another large fraction of all baryons resides within 
mildly-nonlinear structures with temperatures 0.01 — IKeV, called Warm-Hot Intergalactic 
Medium (WHIM). The amount of baryons within this WHIM is estimated to reach 20 - 40% 
(Cen & Ostriker 1999; Dave et al. 1999, 2001). 

Being highly ionized the WHIM could generate temperature anisotropics by means of 
the thermal or kinematic Sunyaev-Zeldovich effect (TSZ, KSZ; Sunyaev & Zeldovich 1972, 
Sunyaev & Zeldovich 1980), adding to the contribution coming from gas in halos. Hallman 
et al (2007) carried out a numerical simulation and found that about a third of the SZ flux 
would be generated by unbound gas. In Atrio-Barandela & Miicket (2006) -hereafter ABM- 
we estimated the CMB anisotropies generated by the WHIM due to the TSZ effect. In 
this letter we compute the imprint on CMB temperature anisotropies due to the peculiar 
motions of gas in filaments and halos. The relative KSZ contribution with respect to the 
TSZ increases with decreasing temperature and the anisotropies generated could be used to 
trace the gas distribution on large scales. Briefly, in Sec. 2 we show the effect of individual 
filaments on the CMB temperature anisotropies. We extend an early calculation of the 
KSZ probability distribution function due to halos by Yoshida, Sheth & Diaferio (2001) by 
including gas not bound in halos and generalize the work of Juszkiewicz, Fisher & Szapudi 
(1998). In Sec. 3 we compute the radiation power spectrum, showing that it reaches its 
maximum at I ~ 400, an angular scale much larger than that of clusters of galaxies; we also 
show that a statistical analysis including the anisotropies generated by the peculiar motions 
in the WHIM improves the fit to WMAP 3yr data. Finally, we present our main conclusions. 



2. The effect of single filament. 

The temperature anisotropy due to the KSZ effect, along a given hue of sight (los), 
is the number density weighted average of the radial component of the peculiar velocity of 
electrons: AT — Tgar J ne{l){ve,r{l) / c)dl — TVr/c. In this expression, Ug is the electron 
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number density, fe,r is the projected peculiar velocity along the line of sight /, (Tt is the 
Thomson cross section and c the speed of light. The last expression follows when the KSZ 
signal is dominated by a single filament of size L and radial component of its peculiar velocity 
Vr- In this case, r = aTU^L is the optical depth along the filament of average electron density 

He- 

To compute the probability distribution that a filament produces a given KSZ anisotropy, 
we assumed that the matter distribution is well described by a log-normal random field. Coles 
& Jones (1990) showed this approximation to be adequate when matter density perturbations 
are weakly non-linear and the velocity field is still in the linear regime. The probability that 
at any given point and at redshift z there is an overdensity ^ = ne{x, z) /ne{z), with ne{z) 
the average electron density, is: = (^cr^v^Svr)""'^ exp [— (log,^ — yU5)^/(2cr|)] where fis and 

as are the mean and rms of the matter density contrast in the linear regime (i.e., fis = 0). 
For filaments of characteristic scale length L, this distribution is also the distribution of the 
optical depth r at each location. The KSZ contribution includes the effect of all baryons 
along the los, independently of whether they are in halos or in the form of WHIM. 

The velocity field is accurately described as Gaussian Random Field with mean and 
rms fivjCfv Its coherence scale is much larger than that of the density field; the motion of 
galaxies, clusters and the IGM filaments is not determined solely by the matter distribution 
in their vicinity, but from a much larger scale. Under this assumption, we can assume that 
density and velocity fields are independent. At each location, the probability distribution 
function (PDF) that an ionized gas cloud would produce a KSZ of a given amplitude is 
simply the product of the distribution of the optical depth (log-normal) and of the velocity 
(gaussian). If X = r and Y = Vr are two continuous random variables with joint probability 
/x,y(a;, y), then the PDF oi U = XY = rvr is: fu{u) = f /x,y(a;, M/a;)|x|^"'^c/x (Rohatgi, 
1976). When both random variables are independent, then fx,Y = fxfv- The resulting PDF 
is symmetric with respect to the y-axis. In Fig. [T]we plot only the region m > 0; the thick 
solid line represents the PDF of the KSZ effect due to the IGM. For comparison, we also 
plot as a thin solid line the gaussian distribution. Both distributions were normalized such 
that the velocity and linear density fields have zero mean and unit variance [a^ = = 1). 
Then, u is expressed in units of {T)av, where (r) is the average optical depth of filaments 
of a size L. Let us remark that the KSZ distribution has much larger tails than the normal 
distribution. The same PDF was computed by Yoshida et al. (2001) for the halo model. In 
their approach, it was unlikely that a given los will contain more than one halo and their 
PDF was dominated by the more abundant and less massive halos, and only in the tails of 
the distribution is the contribution of the more massive halos significant. In ours, the optical 
depth along each line of sight can be much larger since it includes the effect of baryons in 
halos together with the WHIM in low dense regions, and this effect will raise the tails of the 
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PDF. 

More informative is the Distribution Function F{u) = /"^ fu{u')du'. In Fig. [1] we plot 
the quantity P{u > Uq) = 1 — F(uo) (dashed hne) that represents the probabihty that a given 
filament produces a KSZ signal u equal or larger than uq = {T)vr- The shaded regions [0, 1.7] 
and [0, 6.2] correspond to percentiles 0.9 and 0.99. For a gaussian distribution, the probability 
that a random variable has an amplitude 2.7 times the rms deviation is 1%. For the KSZ 
distribution, at the same percentage, the amplitude is 6.2 times the average; if a filament 
of size L and electron density He has an optical depth r ~ 2 x 10~'^(ne/10^m~^)(L/lMpc) 
and moves with a velocity of 300km/ s, it will typically give rise to ATksz ~ 5.5yuK. About 
9% of the filaments in this population will have enough velocity and density to produce an 
effect in the range [9.3, 34]/xK and 1% will have values larger than 34yuK. The same region 
will produce, on average, a thermal SZ signal of 6tsz ~ ll.5G{h')fj,K{Tx/lkeV), where 
~ —2 in the Rayleigh- Jeans regime. 

Observations with the Very Small Array of the Corona Borealis Supercluster (Genova- 
Santos et al. 2005) have detected an extended (~ 30 arcmin) cold spot without significant 
X-ray emission on the R0SAT-R6 map, termed decrement H. Those authors described the 
likelihood that it could be partly generated by an unknown massive cluster in the background 
or by a filament. They did not include the significant KSZ contribution in their study. For 
example, 1% of all the filaments of size ~ 5Mpc, aligned with the los and with a coherent flow 
at similar scale will produce a Stsz + ^ksz ~ (—150, — 200)/iK, enough to explain the central 
temperature of —230 ± 23yuK of decrement H. Further observations, like those to be carry 
out by the PLANCK satellite, at various frequencies and with better angular resolution, are 
required to quantify the KSZ contribution. 



3. Contribution to CMB anisotropies. 

The KSZ effect is independent of redshift; along a los contributions from individual fil- 
aments will average out since their peculiar velocities will have random directions. Since the 
non- linear density and velocity fields evolve with redshift, so does the contribution of the fil- 
ament network. To estimate the temperature anisotropies generated by the peculiar motions 
of filaments and halos we adapted the formalism developed by Choudhury, Padmanabhan & 
Srianand (2001) and ABM. The correlation function of the KSZ signal is 

0(6) = (^ ^Tksz{x)ATksz{x') ^ ^ ^ f'^^'^ r'\ldl\nSx, z)v,Alx, z)n,{lx' , z')v,^r{lx^ 

-'-o -'- o C J o Jo 

(1) 
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where 9 is the angular separation between the two los, x and x'. The double integration is 
carried out from the origin to the highest redshift Zf where the IGM is still fully ionized. 

The main difficulty in making an analytical estimate is to relate density and velocity 
fields at each location since both random fields have very different scale lengths. If we assume 
they are independent, like in the previous section, we can decompose the velocity field into 
two contributions: the random motion of matter with respect to the center of mass Vr{z), 
and the velocity of the center of mass with respect to a comoving observer, also termed Bulk 
Flow Vb'- Ve{lx,z) = Vb{z) + vr{Ix, z). With this decomposition, the integrals containing 
the random component in eq. ([1]) average out and 
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^(^) = 4 / / dldl'VB{R{z))D,{z)VB{R{z'))D,{z'){n,{lx,z)n,{lx',z')) (2) 

^ J o Jo 

where Vb{R{z)) denotes the mean bulk velocity of a sphere with radius R{z), the comoving 
distance from the observer to a filament at redshift z and D^[z) is the velocity linear growth 
factor. 

Assuming a log-normal distribution for matter and baryons, in the non-linear regime, 
the electron number density ne(x, z) in the IGM is: UBi^, z) = uq^z) exp[(5B(x, z) — A'^^{z)/2], 
where x denotes the spatial position at redshift z and \^{z)\ is the proper distance; 5_b(x, z) 
is the linear baryon density contrast, no{z) = Pb(1 + z^/fiB^p and ps, riip are the baryon 
density and proton mass, respectively. Finally, 

A|WH4(X,.))^B^W/|^J|^^|§^, (3) 

where D{z) = D{z, Q^, Qm) is the matter linear growth factor. The spatial average in eq. ([2]) 
is ABM: 

C{e) = 4 / ^ / ^ dldl'VB{R{z))D,{z)VB{R{z'))D,{z')ne{z)ne{z')F{e, z, z'), (4) 

Jo Jo 



with 

F{e,z,z') = -1 +exp 



D{z)D{z') p P{k)edk sin(fc|x-x'|) 



27r2 Jo [I + xl{z)k^][l + xl{z')k^] A;|x-x' 



(5) 



where |x — x'| denotes the proper distance between two filaments at positions ^{z) and 
x'(z') separated by the angle 9, P{k) is the matter power spectrum and x^, is the Jeans 
length at each redshift, that relates the matter power spectrum to that of baryons (Fang 
et al. 1993): Xh{z) = Hq^ [{2^kBTm{z)) / {?>pmpVtrn{^ + z))]^^"^ . In this last expression, Tm 
is the mean temperature of the IGM, 7 ~ (1.2 — 1.6) is the polytropic index, VL^ is the 
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cosmological fraction of matter density and /i is the mean molecular weight of the IGM. The 
helium weight fraction was Y = 0.24. The factor jT^ defines the Jeans length and fixes 
the smallest scale that contributes to the spectrum. We will take •yTm ~ 2 x lO'^K (Schaye 
et al. 2000). The number density of electrons rig in the IGM can be obtained by assuming 
ionization equilibrium between recombination, photo-ionization and coUisional ionization. At 
the conditions valid for the IGM (temperature in the range 10^ — lO^K, and density contrast 
6 < 100) the degree of ionization is very high. Commonly = ens, with 0.9 < e < 1 
depending on the degree of ionization. For future reference, Xb = 0.3 — O.S/i^^Mpc, depending 
on model parameters (see ABM). 

In eq. HI for each 6 and z the function F{6, z, z') is strongly peaked a.t z = z' and is 
different from zero only within a very narrow range G [z — Az,z + Az]. The effective 
peak width 2Az increases and its maximum amplitude decreases with decreasing angle and 
increasing redshift, simplifying the numerical evaluation of eq. HI in our numerical integration 
we took Az = 0.1. Integration around the peak width was carried out with high accuracy, 
using z, z' intervals of (52; = 5 x 10"'^ and angle separations of 30 arcseconds. Further 
integration (with respect to z only) is carried out with a much smoother function and can 
be performed by a relatively rough spacing in redshift. 

In Fig. [2^, we show the radiation power spectrum obtained using the formalism devel- 
oped before. Since some of the parameters in our model are poorly determined, we considered 
the effect on the radiation power spectrum of: (a) the amplitude of the matter power spec- 
trum measured by cxg, (b) the averaged temperature of the IGM, and (c) the maximum 
overdensity ^ for which the log-normal distribution is a good description of the IGM. We 
took as = (0.8; 0.9), = (1.4 x 10^5. x 10^) K and 7 = 1.3 throughout. We restricted 
all integrations to overdensities ^ < 500. The overall shape of the power spectra is very 
similar for a large range of the parameter space: the spectrum peaks around / ~ 200 — 600 
(shaded area) with a maximum amplitude of 10 — 25/iK. This signal is much larger than 
the one considered by Hajian et al (2007), associated to the free electrons in the halo of the 
Galaxy. The amplitude and angular scale depend on Xh or, equivalently, on T^. The Jeans 
length fixes the scale where the power in the baryon density field falls below that of the DM. 
Decreasing x^ decreases the scale where baryons evolve non-linearly and increases the frac- 
tion of baryons sharing the same coherent motion, increasing the KSZ power spectrum and 
shifting the maximum to smaller angular scales. The power also increases with increasing 
either or ^. If we vary parameters within the range ag G (0.7 — 0.9) and ^ G (50 — 500) the 
amplitude of the power spectrum could change up to a factor ^ 5. has the largest effect: 
the amplitude scales as ~ similarly to the Thermal Sunyaev-Zeldovich contribution (see 
ABM). 
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In Fig. [2]d we plot the contributions to the power spectrum from different redshift 
bins {0 < z < 0.15;dz = 0.01). As remarked above, the differential contribution strongly 
decreases with redshift. The signal is dominated by: (1) the bulk motion on scales of 
300/;.~^Mpc and below and (2) by the structures with largest density contrast. Contributions 
from higher redshifts and lower densities are negligible due to the superposition of filaments 
with randomly oriented velocities with decreasing contributions. Like in our earlier results 
for the TSZ effect, the scale of the maximum is fixed by the Jeans length; at z=0.01, xj, ~ 
AOOh^^kpc corresponds to an angular scale / ~ 300, as shown. 

Comparing the spectrum of Fig. [2] with the results of the KSZ effect obtained using 
the halo model (Molnar & Birkinshaw 2000, Cooray & Sheth 2002) the former shows more 
power (factor 3-10, depending on model parameters) and at different angular scale. While 
the halo model only accounts for the gas in halos, our formalism also incorporates the gas in 
low dense regions. Those regions with intermediate density contrast include the largest gas 
fraction at low redshift (see ABM) when the bulk flow velocities are largest, which boost the 
power with respect to the halo model result. 

Let us consider whether the KSZ contribution described above could explain the ex- 
cess of power in the range / ~ 200 — 600 measured by WMAP 3yr data with respect to 
the best fit model (see Hinshaw et al 2007, Fig. 17). We explored the parameter space of 
ACDM models using a Monte Carlo Markov Chain (MCMC); we added a KSZ contribu- 
tion with an amplitude at the maximum of 25/iK to the radiation power spectrum of each 
cosmological model. This approach is not self consistent since when ag changes so does the 
amplitude of the KSZ contribution. Even though, the best fit model improved the fit in one 
unit: Ax^ = 1- This improvement is not statistically significative from a Bayesian point 
of view (Liddle 2004): we introduced two parameters 'yTm and the maximum overdensity 
^, to obtain only a mild improvement on the fit. Comparison with the concordance model, 
parameters deduced from our MCMC differ by less than 1%, being the variation on flfyh'^ 
the most significative. This exercise does show that ionized gas could originate measurable 
CMB temperature anisotropies and CMB observations could provide evidence of the gas 
distribution at large scales. 



4. Discussion. 

If the 'missing baryons' are in the form of a highly ionized WHIM, and this medium 
shares the same velocity field than galaxies and dark matter, then it will produce a contri- 
bution to the CMB temperature anisotropies that could be detectable. Gas in halos and 
filaments with bulk flows aligned with the los could explain a significant fraction of the sig- 
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nal measured in the decrement H of the Corona Boreahs supercluster observed by the VSA 
interferometer. Also, we computed the effect of the overall population of filaments on CMB 
temperature anisotropies. Their contribution has a caracteristic angular scale at / ~ 400, 
very different from that of clusters of galaxies. Uncertainties in our model parameters trans- 
late into an uncertainty of a factor 5 in the amplitude of the power spectrum. 

Intrinsic and KSZ temperature anisotropies have the same frequency dependence, they 
can not be separated directly from the data. Using a MCMC analysis we have shown that 
this extra component improves the fit of the concordance model to WMAP Syr data, but 
the improvement can not be considered significant. WMAP Syr data or the future PLANCK 
observations could have enough statistical power to indicate the existence of this component. 
Further, PLANCK has enough angular resolution and frequency coverage as to measure the 
KSZ contribution of structures such as decrement H in the VSA data of the Corona Borealis 
Supercluster. 
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Ciencia (project FIS2006-05319) and Junta de Castilla y Leon (project SA010C05). We 
thank R. Sheth for discussions. J. P.M. thanks the University of Salamanca for its support 
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Fig. 1. — KSZ (thick solid) and Gaussian (tliin solid) probability distribution functions. The 
dashed line represents P{u > uo), the probability that the a given filament produces a KSZ 
anisotropy u equal or larger than a given value uo- Shaded areas correspond limit the regions 
of these probabilities being 10% and 1%. 
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Multipole I 

Fig. 2. — (a) KSZ Radiation power spectra for different parameter sets cxg, T^, C (see labels), 
(b) Differential Power spectrum contribution from subsequent redshift bins of width dz = 
0.01. The shaded area indicates the interval where power spectra with different parameters 
reach their maxima. 



